function [F,G,H,L,M,N] = pyramst_uhlig(y, x, parms)
G = zeros(2,2);
G(1,1) = 1;
G(1,2) = -parms(1)*parms(5)*y(1)^(-2 + parms(1))*y(2)^(-parms(3))*(1 + y(3))*(1 - parms(1))/(1 + parms(2));
G(2,1) = 1;
G(2,2) = parms(3)*y(2)^(-1 - parms(3));
F = zeros(2,2);
F(1,1) = 0;
F(1,2) = 0;
F(2,1) = 0;
F(2,2) = -parms(3)*y(2)^(-1 - parms(3))*(1 - parms(4) + parms(1)*parms(5)*y(1)^(-1 + parms(1))*(1 + y(3)))/(1 + parms(2));
H = zeros(2,2);
H(1,1) = -1 + parms(4) - parms(1)*parms(5)*y(1)^(-1 + parms(1))*(1 + y(3));
H(1,2) = 0;
H(2,1) = 0;
H(2,2) = 0;
M = zeros(1,2);
M(1,1) = -parms(5)*y(1)^parms(1);
M(1,2) = 0;
L = zeros(1,2);
L(1,1) = 0;
L(1,2) = parms(1)*parms(5)*y(1)^(-1 + parms(1))*y(2)^(-parms(3))/(1 + parms(2));
N = zeros(1,1);
N(1,1) = 0;
